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Abstract 

The accretion of matter onto a massive black hole is believed to feed the relativistic 
plasma jets found in many active galactic nuclei (AGN). Although some AGN accelerate 
particles to energies exceeding 10 12 electron Volts (eV) and are bright sources of very-high- 
energy (VHE) 7-ray emission, it is not yet known where the VHE emission originates. Here 
we report on radio and VHE observations of the radio galaxy M 87, revealing a period of 
extremely strong VHE 7-ray flares accompanied by a strong increase of the radio flux from 
its nucleus. These results imply that charged particles are accelerated to very high energies 
in the immediate vicinity of the black hole. 

Active galactic nuclei (AGN) are extragalactic objects thought to be powered by massive 
black holes in their centres. They can show strong emission from the core, which is often 
dominated by broadband continuum radiation ranging from radio to X-rays and by substantial 
flux variability on different time scales. More than 20 AGN have been established as VHE 7-ray 
emitters with measured energies above 0.1 tera electron Volts (TeV); the jets of most of these 
sources are believed to be aligned with the line-of-sight to within a few degrees. The size of the 
VHE 7-ray emission region can generally be constrained by the time scale of the observed flux 
variability El 13 but its location remains unknown. 

We studied the inner structure of the jet of the giant radio galaxy M 87, a known VHE 7-ray 
emitting AGN EH ED 13 with a (6.0 ± 0.5) x 10 9 M Q black hole HI, scaled by distance, located 
16.7 Mpc (54 million light years) away in the Virgo cluster of galaxies. The angle between its 
plasma jet and the line-of-sight is estimated to lie between 15 — 25 deg (see supporting online 
text). The substructures of the jet, which are expected to scale with the Schwarzschild radius 
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R s of the black holcQ are resolved in the X-ray, optical and radio wavebands (Fig. [I]). High- 
frequency radio very long baseline interferometry (VLBI) observations with sub-milliarc second 
(mas) resolution are starting to probe the collimation region of the jet [81 . With its proximity, 
bright and well-resolved jet, and very massive black hole, M 87 provides a unique laboratory 
in which to study relativistic jet physics in connection with the mechanisms of VHE 7-ray 
emission in AGN. 

VLBI observations of the M 87 inner jet show a well resolved, edge-brightened structure 
extending to within 0.5 mas (0.04 pc or 70 R s ) of the core. Closer to the core, the jet has a 
wide opening angle suggesting that this is the collimation region [8|. Generally, the core can 
be offset from the actual location of the black hole by an unknown amount [9J, in which case it 
could mark the location of a shock structure or the region where the jet becomes optically thin. 
However, in the case of M 87 a weak structure is seen on the opposite side of the core from the 
main jet, which may be the counter-jet, based on its morphology and length [fTOlfTTTl . Together 
with the observed pattern in opening angles, this suggests that the black hole of M 87 is located 
within the central resolution element of the VLBI images, at most a few tens of R s from the radio 
core (see supporting online text). Along the jet, previous monitoring observations show both 
near-stationary components [11] (pc-scale) and features that move at apparent superluminal 
speeds lfT2"l[l"3l (100 pc-scale). The presence of superluminal motions and the strong asymmetry 
of the jet brightness indicate that the jet flow is relativistic. The near- stationary components 
could be related to shocks or instabilities, that can be either stationary, for example if they are 
the result of interaction with the external medium, or slowly moving if they are the result of 
instabilities in the flow. 

A first indication of VHE 7-ray emission from M 87 was reported by the High Energy 
Gamma-Ray Astronomy (HEGRA) collaboration in 1998/99 [3]. The emission was confirmed 
by the High Energy Stereoscopic System (H.E.S.S.) in 2003-2006 [2], with 7-ray flux variability 
on time scales of days. M 87 was detected again with the Very Energetic Radiation Imaging 
Telescope Array System (VERITAS) in 2007 [4] and, recently, the short-term variability was 
confirmed with the Major Atmospheric Gamma-Ray Imaging Cherenkov (MAGIC) telescope 
during a strong VHE 7-ray outburst [5] in February 2008. Causality arguments imply that the 
emission region should have a spatial extent of less than « 5SR S , where 5 is the relativistic 
Doppler factor. This rules out explanations for the VHE 7-ray emission on the basis of (i) dark 
matter annihilation lfl4ll . (ii) cosmic-ray interactions with the matter in M87 lfl5l . or (iii) the 
knots in the plasma jet (Fig. [Tp). Leptonic [fT6l[T7ll and hadronic [18J VHE 7-ray jet emission 
models have been proposed. However, the location of the emission region is still unknown. 
The nucleus ||T9~1 1201 , the inner jet [1211 or larger structures in the jet, such as the knot HST-1 
(Fig. [lp), have been discussed as possible sites [13J. Because the angular resolution of VHE 
experiments is of the order of 0.1 deg, the key to identifying the location of the VHE 7-ray 
emission lies in connecting it to measurements at other wavebands with considerably higher 
spatial resolutions. An angular resolution more than six orders of magnitude better (less than 

'The Schwarzschild radius of a black hole with the mass m is defined as R s = 2Gm/c 2 , G is the gravitational 
constant, and c is the speed of light. The Schwarzschild radius defines the event horizon of the black hole. 
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6x10 degrees, corresponding to approximately 30 R s in case of M 87) can be achieved with 
radio observations (Fig.[T]). 

We used the H.E.S.S. [22], MAGIC (23 and VERITAS O instruments to observe M 87 
during 50 nights between January and May 2008, accumulating over 95 h of data (corrected 
for the detector dead times) in the energy range between 0.1 TeV and several 10's of TeV. 
Simultaneously, we monitored M 87 with the Very Long Baseline Array ll25l (VLB A) at 43 GHz 
with a resolution of 0.21 x 0.43 mas li26ll . corresponding to about 30 x 60 R s , see ll27l . During 
the first half of 2008, three X-ray pointings were performed with the Chandra satellite [28]. Our 
light curves are shown in Fig. [2} 

We detected multiple flares at VHE in February 2008 with denser sampling, following a 
trigger sent by MAGIC [~23 h of the data published in [5]]. The short-term VHE variability, 
first observed in 2005 [2], is clearly confirmed and the flux reached the highest level observed 
so far from M 87, amounting to more than 10% of that of the Crab Nebula. At X-ray frequencies 
the innermost knot in the jet (HST-1) is found in a low state, whereas in mid February 2008 the 
nucleus was found in its highest X-ray flux state since 2000 11281 . This is in contrast to the 2005 
VHE 7-ray flares [2], which happened after an increase of the X-ray flux of HST-1 over several 
years [|29|| , allowing speculation that HST-1 might be the source of the VHE 7-ray emission 
[TT3l : no 43 GHz radio observations were obtained at that time. Given its low X-ray flux in 
2008, HST-1 is an unlikely site of the 2008 VHE flaring activity. 

Over at least the following two months, until the VLBA monitoring project ended, the 
43 GHz radio flux density from the region within 1.2 mas of the core rose by 30% as compared 
with its level at the time of the start of the VHE flare and by 57% as compared with the average 
level in 2007 (Fig. [2]>. The resolution of the 43 GHz images corresponds to 30 x 60i? s and 
the initial radio flux density increase was located in the unresolved core. The region around the 
core brightened as the flare progressed (Fig. [3]), suggesting that new components were emerging 
from the core. At the end of the observations, the brightened region extended about 0.77 mas 
from the peak of the core, implying an average apparent velocity of 1.1 c (c is the speed of 
light), well under the approximately 2.3 c seen just beyond that distance in the first half of 2007. 
Astrometric results obtained as part of the VLBA monitoring program show that the position 
of the M 87 radio peak, relative to M 84, did not move by more than ~6 R s during the flare, 
suggesting that the peak emission corresponds to the nucleus of M 87. 

Because VHE, X-ray and radio flares of the observed magnitude are uncommon, the fact 
that they happen together (chance probability of P < 0.5%, supporting online text) is good 
evidence that they are connected. This is supported by our joint modeling of the VHE and radio 
light curves: The observed pattern can be explained by an event in the central region causing 
the VHE flare. The plasma travels down the jet and the effect of synchrotron self-absorption 
causes a delay of the observed peak in radio emission because the region is not transparent 
at radio energies at the beginning of the injection (supporting online text, Sec. 3). The VLBI 
structure of the flare along with the timing of the VHE activity, imply that the VHE emission 
occurred in a region that is small when compared with the VLBA resolution. Unless a source 
of infrared radiation is located very close to the central black hole, which is not supported by 
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current observations 11301 . TeV 7-ray photons can escape the central region of M87 without 
being heavily absorbed through e + e~ pair production [fT9ll20l . 

The light curve might indicate a rise in radio flux above the range of variations observed 
in the past, starting before the first VHE flare was detected. This could imply that the radio 
emission is coming from portions of the jet launched from further out in the accretion disk than 
that responsible for the VHE emission. However, it is difficult to derive a quantitive statement 
on this, because no VHE data were taken in the week previous to the flaring. Thus, an earlier 
start of the VHE activity cannot be excluded, either. 

A possible injection of plasma at the base of the jet observed at optical and X-ray energies 
with a delayed passage through the radio core ~10 4 R s further down the jet - interpreted as 
a standing shock and accompanied by an increase in radio emission - has been discussed in 
the case of BLLac flU (with evidence for VHE emission, see supporting online text for more 
details). M 87 is much closer than BLLac and has a much more massive black hole, allowing 
the VLBA to start resolving the jet collimation region whose size, from general relativistic 
magnetohydrodynamic simulations [|3TI . is thought to extend over ~1000 R s . In case of M 87 
the radio core does not appear to be offset by more than the VLBA resolution of ~50 R s from 
the black hole (see supporting online text) and the jet has a larger angle to the line-of-sight than 
in BLLac. Thus the coincidence of the VHE and radio flares (separated in photon frequency by 
16 orders of magnitude), constrains the VHE emission to occur well within the jet collimation 
region. 
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Figure 1 : M 87 at different photon frequencies and length scales. (A): Comparison of the different length scales. 
(B): 90 cm radio emission measured with the VLA. The jet outflows terminate in a halo which has a diameter 
of roughly 80kpc (15')- The radio emission in the central region is saturated in this image. Credit: F.N. Owen, 
J. A. Eilek and N.E. Kassim ||32ll , NRAO/AUI/NSF. (C): Zoomed image of the plasma jet with an extension of 
2 kpc (20"), seen in different frequency bands: X-rays (Chandra, upper panel) optical (V band, middle) and radio 
(6 cm, lower panel). Individual knots in the jet and the nucleus can be seen in all three frequency bands. The 
innermost knot HST-1 is located at a projected distance of 0.86 arcseconds (60 pc, w 10 5 R s ) from the nucleus. 
Credit: X-ray: NASA/CXC/MIT/H. Marshall et al., radio: F. Zhou, F. Owen (NRAO), J. Biretta (STScI), optical: 
NASA/STScI/UMBC/E. Perlman et al., |7|. (D): An averaged, and hence smoothed, radio image based on 23 
images from the VLBA monitoring project at 43 GHz. The color scale gives the logarithm of the flux density in 
units of 0.01 mJy/beam. The indication of a counter-jet can be seen, emerging from the core towards the lower 
left side. Image from ll27l . 
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Figure 2: Combined M 87 light curves from 2007 to 2008. (A): VHE 7 -ray fluxes (E > 0.35 TeV, nightly av- 
erage), showing the H.E.S.S., MAGIC and VERITAS data. The fluxes with statistical errors (1 standard deviation) 
were calculated assuming a power-law spectral shape of dN/dE oc E~ 23 . Monthly-binned archival VERITAS 
data taken in 2007 are also shown |4|. The systematic uncertainty in the flux calibration between the experiments 
was estimated to be on the order of 20% based on Crab Nebula data. The regular gaps in the light curve correspond 
to phases of full moon during which no observations were possible. The inlay shows a zoomed version of the flar- 
ing activity in February 2008; the time span is indicated by the grey vertical box in all panels. (B): Chandra X-ray 
measurements (2 — 10 keV) of the nucleus and the knot HST-1 [28 1. (C): Flux densities from the 43 GHz VLBA 
observations are shown for (i) the nucleus (circular region with radius r = 1 .2 mas = 170 R s centered on the peak 
flux), (ii) the peak flux (VLBA resolution element), and (iii) the flux integrated along the jet between distances of 
r = 1.2 — 5.3 mas (compare with Fig. [3]). The error bars correspond to 5% of the flux. The shaded horizontal area 
indicates the range of fluxes from the nucleus before the 2008 flare. Whereas the flux of the outer regions of the 
jet does not change substantially, most of the flux increase results from the region around the nucleus. Image from 
ED. r 
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Figure 3: VLBA images of M87 at 43 GHz. (A): Average (hence smoothed) of 11 images from data taken 
between January and August, 2007, well before the VHE and radio flare. The contour levels start at 5, 10, 14.3, 
and 20 mJy per beam and increase from there by factors of y/2. The restoring beam used for all of the images 
is 0.21 x 0.43 mas (30 x 60i? s ) elongated in position angle — 16deg, as shown by the ellipse in the upper left 
corner. (B): Image from 5 April 2008, with the same contours and colors as in (A). The linear scale in light days 
and Schwarzschild radii is also shown. (C - F): Difference images for observations during the period of the radio 
flare showing its effects. These were made by subtracting the average image (A) from the individual epoch images. 
The contours are linear with 10 (white) at intervals of 7mJy per beam followed by the rest (black) at intervals of 
70 mJy per beam; negative contours are indicated by dashed lines. At the time of the VHE flare, the core flux 
density was already above the average but the region of the jet between —0.5 and — l.Omas RA offset was below 
average, suggesting that there had been a period of below-normal activity leading up to the flare and that the radio 
flare may have begun before the VHE flare. The sequence shows the substantial rise in the core flux density and 
the appearance of enhanced emission along the inner jet. Image from J27). 
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Supporting Online Material 



1) The geometry of the M 87 jet 

General assumptions. The luminosity distance of M 87 is D = 16.7 Mpc [|33l, so that an 
angle of 1 milli-arc sec (mas) corresponds to 0.081 pc = 2.5 x 10 17 cm. The mass of the central 
black hole has been recently determined using detailed modeling of long slit spectra of the 
central regions of M 87 to be (6.0 ± 0.5) x 10 9 M BH, corrected for the distance we use^J The 
Schwarzschild radius is i? s ~ 1.8 x 10 15 cm; lmas therefore corresponds to approximately 
140 R s . The observed timing delay At obs between the onset of the VHE flare and the radio peak 
is about 50 days (~4x 10 6 s). The following equations hold: The observed velocity in units of 
the speed of light c is calculated as /3 obs = [/3 int sin 9]/[l — (3- mt cos 9], the Doppler boost factor 
is 5 = [r(l — /3i nt cos^)] -1 , the observed time delay is At obs = A£ int [l — (3- mt cos^)], and the 
distance along the jet is r int = r obs / sin 6 ps 3 r obs for 9 = 20 deg (see below). 

The jet orientation angle. The most rapid optical proper motions are observed along the jet 
with superluminal speeds lfi~2l of /3 obs ~ 6 (HST-1), /3 obs ~ 5 (knot D), and /3 obs ~ 4 (knot E). 
On the other hand, the most rapid radio proper motions seen from features moving through HST- 
1 indicate superluminal speeds in the range of /3 obs ^ 3 — 4 [13] and through knot D indicate 
a superluminal speed of (3 ohs m 2.5 11331 . Assuming the jet plasma moves at the speed of light 
(/? int = 1), one obtains a maximum viewing angle of cos(6 l max ) = [(3% hs — l]/[/3^ hs + 1]. The 
highest optically determined superluminal speed requires a jet viewing angle of 9 < 18.9 deg. 
On the other hand, the highest radio determined superluminal speed requires 9 < 28.0—36.9 deg 
for (3 ohs w 4 — 3, respectively. Jet angles of 9 = 30 — 45 deg were derived IflOll based on radio 
observations at 43 GHz, assuming the component to the East of the core is a counter-jet, and 
with a velocity measurement based on only one pair of observations. In this paper we assume a 
likely range of the jet angle of 9 = 15 — 25 deg. 

The jet opening angle. The observations indicate the following jet full opening angles as a 
function of the distance to the core: VWs(< 0.5 mas) ps 60 deg flU, VWs(l — 5 mas) ps 12 deg 
[|26ll , and ^obs(> 10 mas) ps 6 deg IfTTTl . Assuming a jet viewing angle of 9 = 20 deg the 
corresponding intrinsic full opening angles are: VVt(< 0.5 mas) ~ 30 deg, ip- mt (l — 5 mas) ~ 
4 deg, and V'mt(> 10 mas) ~ 2 deg. The observed pattern of opening angles suggest that the 
radio core corresponds to the formation region of the jet. 

The counter-jet. The optical identification of an emission feature observed at radio frequen- 
cies located 24 arc sec away from the nucleus in the direction opposite to the jet resulted in a 

2 Gebhardt & Thomas (2009) assumed a distance of 17.9 Mpc. As those authors point out, that mass is about 
70% higher than previous determinations J34), corrected for the different assumed distances. The reason for such 
a large change is not fully explained yet. The use of the new mass |6] implies a higher resolution of our VLBA 
observations in terms of R s , but does not substantially affect our conclusions. 
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first indication of a counter-jet in M87 ll36ll . This feature is also seen in observations at mid- 
infrared frequencies [ 30 1 . Radio observations at wavelengths of 2 cm show clear indications of 
a counter-feature extending up to 3 mas from the radio core in the direction opposite to the jet 
ffTTTl . While individual features move along the jet (determined from 21 images taken over a time 
span of more than 10 years), the counter feature seems to move in the opposite direction with 
an apparent velocity of (0.010 ±0.001) c, strengthening the counter-jet interpretation. However, 
it cannot be fully excluded that this apparent movement is a result of temporal under- sampling 
in the kinematic analysis, although the individual tracked jet features give consistent bright- 
ness temperatures across the observation epochs, which would be unlikely in case of incorrect 
cross-identifications of the components. The jet to counter-jet brightness ratio (0.5 — 3.1 mas) 
is 10 — 15. A counter-feature is also observed at 22 and 43 GHz, extending roughly 1 mas from 
the core in the direction opposite to the jet ifTOl . see also Fig. [lp. The jet to counter-jet bright- 
ness ratio is calculated to be ~ 14. Marginal indication for an apparent velocity of 0.17 c of 
the counter-jet away from the core is found from three 43 GHz radio images. The observations, 
however, suggest a temporal under- sampling of those data. 

Position of the Black hole. At radio frequencies of 86 GHz, the core is no larger than 25 x 7 R s 
113711 . and its jet morphology is consistent with a wide-opening angle jet base, as seen at 43 GHz, 
converging near the point of the black hole. Farther out between about 1 and 10 mas the bright 
edges of the jet converge toward a point close to the maximum extent of the counter-feature. 
Interpreting it as a part of the jet would require that the jet before the radio peak is better 
collimated than afterwards, although it appears to be as wide as the jet itself. The astrometric 
results relative to M 84 show that the root mean square scatter of the position of the radio core 
is about 5 x 2 R s along and across the beam with no clear systematic motion (Davies et al., in 
preparation). If the radio core is farther down the jet and the jet power is going up fractionally 
like the flux density, this would imply a very stable position of the shock region. The spatial 
stability is a reasonable assumption for the radio core being located close to the black hole. 
Theoretical modelling also supports the hypothesis that the 43 GHz radio peak emission results 
from the position of the black hole ll38l . 

2) The VHE and radio observations and frequency of flares 

VHE. The H.E.S.S., MAGIC and VERITAS collaborations operate imaging atmospheric Che- 
renkov telescopes (IACTs) located in Namibia, the Canary Islands (Spain) and Arizona (USA), 
respectively. The telescopes measure cosmic 7-ray photons (entering the atmosphere of the 
Earth) in an energy range of 0. 1 TeV up to several 10's of TeV. M 87 has been observed at those 
energies for the last ten years. Except for the 2008 observation campaign, the observations were 
scheduled in advance and did not follow any external or internal triggers, leading to arbitrarily 
sampled light curves. During the observations of the past 10 years only two episodes of flaring 
activity have been measured: in 2005 [2] and in 2008 (reported in this paper). For the first time, 
M 87 was observed by H.E.S.S., MAGIC and VERITAS in a joint campaign for more than 120 h 
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in 2008 (more than 95 h of data after quality selection). The integral fluxes presented in this 
paper (Fig. [2} were calculated^ under the assumption that the spectrum of M 87 is described by 
a power-law function dN/dE oc E -2 3 Q. Any correlation between the spectral shape and the 
flux level has not yet been established for M 87. The relative frequency of flaring activity was 
estimated by fitting the night-by-night binned light curves as measured by H.E.S.S., MAGIC 
and VERITAS with a constant function (using all available data - partly archival - from 2004 
to 2008). Subsequently flux nights with the most significant deviation from the average were 
removed until the fit resulted in a reduced x 2 P er degree of freedom of less than 1; all removed 
points corresponded to flux values higher than the average. The light curves are compatible with 
constant emission for 49 out of 53 nights (H.E.S.S., 2004-2008), 12 out of 21 nights (MAGIC, 
2008) and 50 out of 51 nights (VERITAS, 2007-2008). Combining these numbers one finds 
flaring activity in the so far recorded data in 14 out of 125 nights of observations, resulting in 
a relative frequency of flares on the order of 10% of all observed nights. Almost all data were 
recorded arbitrarily and except for four nights (with a time difference of ~ 0.5 days between 
the VERITAS and the H.E.S.S./MAGIC observations) all observations were separated in time 
by more than one day. Therefore we assume that this number gives an estimate of the general 
chance to measuring a VHE 7-ray flare from M 87. However, the relative frequency of flaring 
activity is overestimated by the fact that the 2008 observations were intensified for some nights 
during the high flux state following the VHE trigger by MAGIC ||5). 

X-ray. M87 was regularly observed at X-ray energies with Chandra, resulting in 61 mea- 
surements of the X-ray flux of the nucleus during the last ten years ll28l . Three measurements 
exceed a flux level higher than 2 times the root mean square (RMS) of the average flux of all 
data points (relative occurrence of ~ 5%). Only one measurement exceeds the level of 3 RMS 
which was taken during the radio flare with a deviation of ~ 4.3 RMS (Fig.[2]in the main text). 

Radio. Throughout 2007, M 87 was observed with the VLBA on a regular basis roughly every 
three weeks If39~l . The aim of this 'movie project' was to study morphological changes of the 
plasma jet with time. Preliminary analysis of the first 7 months showed a fast evolving structure, 
somewhat reminiscent of a smoke plume, with apparent velocities of about twice the speed of 
light. These motions were faster than expected so the movie project was extended from January 
to April 2008 with a sampling interval of 5 days. A full analysis of these data is in progress 
and details will be published elsewhere. The observed radio flux densities reached at the end 
of the 2008 observations, roughly 2 months after the VHE flare occurred, are larger than seen 
in any previous VLBI observations of M 87 at this frequency, including during the preceding 
12 months of intensive monitoring, in 6 observations in 2006 and in individual observations in 
1999, 2000, 2001, 2002, and 2004 El. The MOJAVE project web siteQgives 15 GHz VLBA 

3 the H.E.S.S. flux points - measured with a higher energy threshold - were extrapolated down from ~ 1 TeV. 
4 http://www. physics.purdue.edu/MOJAVE/ 
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flux densities at 27 epochs since 1995, with the highest flux value measured on May 1, 2008^J 
most of the data (except the last two epochs) are published in [|40i . Assuming a flare duration 
of ~4 months, a similar flare was not observed during a total period of 5 x 4 = 20 months 
based on the 5 observations from 2004 and earlier, which are well separated, 8 months taking 
into account overlap based on the 2006 pilot observations spread over 4 months, and 14 months 
during the 2007/2008 monitoring including 2 months before the start but not including the time 
during the observed flare. That is 42 months total for which a similar flare was not in progress. 
By the same accounting, there are 4 months with a flare. So the probability of a radio flare 
being in progress at any given time is 4/46 ~ 10%, suggesting that radio flares of the observed 
magnitude are uncommon. 

The observed radio/X-ray/VHE-pattern. The probability of observing k out of n nights ex- 
ceeding a flux baseline (for which the chance probability is p) is: p(n, k) = On ■ p k ■ (1 —p) n ~ k . 
The probability of observing k or more flare nights is P(n, > k) = Y,]=kP( n i j)- A joint esti- 
mated of the chance probability of the observed pattern is difficult, since the characteristic time 
scales of flux changes are different and the data are not sampled equally, or are partly sam- 
pled with higher frequencies as the characteristic time scale of the flux changes (over- sampling) 
or much less frequent (under- sampling). However, if defining a time window covering the 
whole increase of the radio flux between end of January to mid of April one finds 8 out of 40 
VHE measurements (Pvhe) and 1 out of 2 X-ray measurements (Px-ra y ) exceeded their base- 
lines. The chance probability of observing this pattern during the radio flare is on the order of 
P = p VHE . P x _ ray = 0.026 • 0.095 < 0.5%. The fluxes in the radio, X-ray, and VHE bands 
reached their highest archival level during the defined time window. 

3) Time-dependent modelling of the radio emission as synchrotron emission 
from a slow outer sheath 

The model. Whereas the VHE 7-ray emission from M 87 varies on time scales of a few days, 
the 43 GHz radio emission from the nucleus steadily increases over a time period of two months 
(Fig. [2]>. A model calculation was performed to test if the slow variations of the radio flux can 
be explained with a self-absorbed synchrotron model of electrons injected into a "slow outer 
sheath" of jet plasma. The slow outer sheath has the geometry of a hollow cone, an assumption 
which is supported by the edge-brightened structure of the jet observed at radio frequencies 
(Fig. [Tp)). As the plasma travels down the jet, it expands, leading to a decline of the frozen-in 
magnetic field and to adiabatic cooling of the electrons. In the model, a 7-ray flare leads to the 
injection of radio-emitting plasma at the base of the jet and at the base of the slow outer sheath. 
The model assumes that the VHE 7-rays are produced very close to the black hole, e.g. in the 
black hole magnetosphere and does not attempt to describe the production mechanism. It merely 

5 This is only a few weeks later than the maximum of the 43 GHz flux reported in this paper. However, no 
further 15 GHz data for 2008 are listed on the MOJAVE web page. 
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Figure 4: Simulated 43 GHz light 
curve resulting from a single injection 
of radio-emitting plasma at time t = 
<T = 1.01, /3 jc t = 0.14, B = 0.5 G). 
Initially, the plasma is opaque owing 
to synchrotron self-absorption. As the 
plasma travels away from the point of 
injection, it expands and becomes trans- 
parent. The expansion leads to a de- 
crease of the magnetic field and to adi- 
abatic cooling of the electrons, and thus 
to the decline of the radio emission. 
The small deviations during the injec- 
tion stage are due to numerical effects. 



uses the observed 7-ray fluxes to normalize the energy spectrum of the electrons responsible 
for the radio emission. In the beginning, the radio-emitting plasma is optically thick, and the 
synchrotron emission cannot escape. Owing to the adiabatic expansion, the plasma eventually 
becomes optically thin leading to a radio flare. The radio flare dies down owing to the decline 
of the magnetic field and the adiabatic cooling of the electrons. 

Following the injection of radio-emitting plasma at time t , a ring of plasma with radius 
R = R s + Pj e tc sin a(t — t ) (with the thickness of the radio bright sheath being l/5th of 
the cone radius) travels down the jet. The emission of the ring is computed in the frame of 
the moving plasma, assuming that the magnetic field scales as B oc 1/R, and the electrons 
cool adiabatically. The calculation uses the standard equations for Lorentz transforming the 
emission of different sections of the ring, see for example IflTTl . Taking into account light travel 
time effects, the received radio flux at 43 GHz is computed. The overall normalization of the 
electron energy spectrum is adjusted to reproduce the observed radio flux. The results of the 
model strongly depend on the choice of the minimal radius -R m i n = -R s , and weakly depend on 
the choices of the magnetic field B, the jet opening angle a, and and the thickness of the sheath. 

Results. We assume an intrinsic cone opening angle of a = 5 deg, a jet angle of 9 = 20 deg 
and a rather low magnetic field of B = 0.5 G at the base of the jet for which radiative cooling can 
be neglected. The simulated radio light curve produced by a single injection of radio-emitting 
plasma is shown in Fig.[4]for an assumed bulk Lorentz factor of T = 1.01 (/3j et = 0.14, giving 
the best fit result). The radio flux needs approximately 20 days to reach its maximum. Figure [5] 
shows the corresponding radio light curve obtained when choosing a time-dependent electron 
injection function proportional to the measured VHE 7-ray fluxes, starting with the VHE data 
taken in January 2008. The spatial extent of the predicted radio source after 50 days is ~ 3 light 
days and therefore still within the central resolution element of the VLBA observations (Fig. [3]). 
The model was chosen to minimize the number of free parameters and assumptions. Other de- 
pendencies could affect the results as follows: (i) If the emitting plasma is more compact (and 
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Figure 5: Observed (data points) 
and modelled (line) 43 GHz radio flux. 
The modelled curve was obtained by us- 
ing the measured VHE 7-ray light curve 
as source function for injecting radio- 
emitting plasma into the slow sheath of 
the jet. The Bulk Lorentz factor used 
is T = 1.01 G9jet = 0.14). No VHE 
7-ray observations were possible during 
phases of full moon, leaving some un- 
[mjd] certainties for the injection function. 

the volume filling factor is < 1), it stays synchrotron self-absorbed for a longer time, increasing 
the time lag between the 7-ray flares and the rise of the radio flux, (ii) The emission volume may 
expand slower than proportional to t 2 , as assumed in the model. Slower expansion would slow 
down the time scale for the rise and decay of the radio flux, (iii) A higher value of /3j et would 
result in a faster radio flare, (iv) All non-thermal particles are injected into the jet right at the 
base of the slow sheath. This model assumption may not be accurate. Additional non-thermal 
particles may be accelerated further downstream which would lead to a longer duration of the 
radio flare, (v) If the magnetic field at the base of the jet is stronger, radiative cooling is not neg- 
ligible any longer and fitting the data would require an assumption of continued acceleration as 
the outer shell flows down the jet. (vi) In a turbulent jet flow, efficient stochastic re-acceleration 
may occur, which could change the picture. However, the current model calculations show that 
synchrotron self-absorption may play a role in explaining the observed slower turn-on of the 
radio emission. 

4) Investigation of VHE 7-ray models for the M 87 VHE/radio flare 

The key questions for the understanding of the VHE 7-ray emission measured from the radio 
galaxy M 87 and from the more than 20 known VHE 7-ray blazars are: (i) What is the under- 
lying particle distribution which is accelerated, (ii) what are the mechanisms to generate the 
7-rays, and (iii) where is the region of the emission located. In the following paragraphs a se- 
lection of models discussed for M 87 in the literature is investigated with a focus on the question 
whether they can explain the observed VHE/radio light curves. Note, however, that some of the 
models have difficulties in explaining the observed hard VHE 7-ray spectra [EMU, which will 
however not be discussed in more detail. 

4.1) Black hole magnetosphere models 

Models of VHE emission from the black hole magnetosphere. The electromagnetic mech- 
anism of extraction of the rotational energy of black holes by the Blandford & Znajek scenario 
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[|42ll seems to be a viable mechanism for powering the relativistic jets of AGN. Particles can 
be accelerated by the electric field of vacuum gaps in the black hole magnetosphere |[T9ll (the 
electric field component parallel to the ordered magnetic field is not screened out) or due to 
centrifugal acceleration in an active plasma-rich environment, where the parallel electric field is 
screened lEOll . Synchrotron and curvature radiation of the charged particles, and inverse Comp- 
ton scattering of thermal photons can produce VHE 7-ray photons 11431 . An important question 
in this scenario is if the 7-rays can escape the central region or if they are absorbed through pair 
creation processes with either photons from the accretion disk IfTBl or infrared photons emitted 
by a potential dust torus for which no clear observational evidence is found so far [|30ll44l . If 
M 87 harbours a non-standard (advection-dominated) accretion disk, 7-rays could escape with- 
out being absorbed ll20ll . An alternative scenario could be that the primary photons create a pair 
cascade whose leakage produces the observed 7-ray emission [45]. The delayed radio emission 
could be explained by the effect of synchrotron self-absorption (see Sec. 3) or the time needed 
to cool the electrons before they dominantly emit synchrotron radiation in the radio regime. 

4.2) Hadronic jet models 

The model of Reimer et al. (2004). In this model a primary relativistic electron population is 
injected together with high-energy protons into a highly magnetized emission region [18]. The 
VHE emission is dominated by either / u ± /7r ± synchrotron radiation or by proton synchrotron 
radiation. The low-energy component is explained by the synchrotron emission of the electron 
population. However, the radio flux is underestimated by the (steady-state) model (explanations 
are discussed in the paper) so that a discussion of the observed radio/VHE flare is beyond the 
scope of this particular model. 

4.3) Leptonic jet models 

The model of Lenain et al. (2008). In this model the high-energy emission region (X-rays up 
to VHE) consists of small blobs (~ 10 14 cm) travelling through the extended jet and radiating at 
distances just beyond the Alfven surface ifTTl . The emission takes place in the broadened jet for- 
mation region, in the innermost part of the jet (corresponding to the central resolution element 
in the 43 GHz radio map). This multi-blob model is a two-flow model, where the fast, compact 
blobs contribute to X-rays and 7-rays through the synchrotron self-Compton mechanism, and 
are embedded in an extended, diluted and slower jet emitting synchrotron radiation from radio 
to optical frequencies. Even though this model describes only steady state emission, the ob- 
served radio/VHE variability can be discussed qualitatively. For instance, a sudden rise of the 
density of the underlying leptonic population at the stationary shock (i.e. in the blobs) translates 
into a flare of X-rays and 7-rays, but no immediate rise of radio emission is expected, because 
the emission volume is synchrotron self-absorbed at radio frequencies, see Sec. 3. However, as 
the flare propagates into the extended, less magnetized, neighbouring jet, the leptons in the jet 
are energized and could cool by emitting at radio frequencies with some delay, creating a di- 
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luted radio flare in response to the VHE flare. An alternative scenario by Giannios et al. (2009) 
explains the fast variability of VHE 7-ray radiation in blazars as a possible result from large 
Lorentz factor (100) filaments within a more slowly moving jet flow [|46ll . 

The model of Tavecchio and Ghisellini (2008). According to this model [12111 the jet consists 
of a fast spine and slow sheath layer. The photons from the fast spine are external-Compton 
boosted to VHE by the slower sheath. In this framework one may assume that the VHE flare 
comes from near to the core and the time lag to the radio maximum is entirely a result of 
propagation of some disturbance down the jet and the associated reduction in synchrotron self- 
absorption. The X-rays should be synchrotron emission that is not self-absorbed and flare at 
about the same time as the VHE 7-rays. In the radio data we see significant edge brightening 
suggestive of a de-boosted spine. Note, however, that edge brightening like that observed can 
also be produced by enhanced surface emissivity. In the model the radio emission originates 
from a region different from that producing the VHE emission, so that a strict flux correlation 
is not required. Detailed modeling would be needed to explain the observed light curves in the 
framework of this model. 

The model of Georganopoulos et al. (2005). In this model the jet decelerates over a length 
of 0.1 pc = 3 • 10 17 cm [16]. The VHE emission is assumed to come from the fast moving 
part near the jet base by inverse-Compton scattering of low-frequency photons from the slower 
moving part of the jet. The model calculations are steady state so that they are difficult to apply 
to an ejection event. However, one can assume that the VHE flare is directly associated with the 
ejection event. As the disturbance propagates down the jet and decelerates we see the radio rise 
later as a result of synchrotron self-absorption effects, similar to the model described in Sec. 3. 
Here the X-ray flux can still rise and fall with the VHE if it comes directly from the disturbance; 
the majority of the observed power comes from the slower part of the flow, but this is a more or 
less steady state result. A similar model by Levinson (2007) of radiative decelerating blobs in 
the jets of VHE 7-ray emitting blazars is described in 071. 

4.4) Jet base / standing shock models 

The model of Marscher et al. (2008). The model is based on the observation of a double 
flare in BLLac BUl. The first flare is seen at X-ray /optical energies accompanied by polariza- 
tion measurements at a date of 2005.82. It indicates an injection event into the jet acceleration 
and collimation region near to the black hole. The flare was followed by the appearance of 
a new radio component in VLBA images that approached and passed through the radio core, 
accompanied by a second X-ray /optical flare at 2005.92, where the 14.5 GHz radio flux begins 
to increase and peaks at about 2006.0. This is interpreted as the passage of the disturbance 
through a standing shock in the jet at a distance of ~ 10 4 R s , see Fig. 3 in BH. VHE 7-ray emis- 
sion has been detected in 2005, however, no evidence for flux variability has been found in the 
data. In this picture the radio core is located at the standing shock and the radio emission in the 
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acceleration and collimation region may be (i) either intrinsically weak or (ii) synchrotron self- 
absorbed. The non-coincidence between the radio peak and the second peak in the X-ray/optical 
is proposed to arise from the longer lifetime of particles radiating at radio frequencies. Since 
the disturbance passes down the expanding jet the radio emission might last longer than the 
second X-ray/optical flare but does not increase in strength after the disturbance passes through 
the shock. 

Applying this model to M 87 one can assume that the observed VHE flare corresponds to 
(A) the first or (B) the second flare. A: The VHE flare indicates the injection at the base of the 
jet and the increase of the radio flux corresponds to the passage through a standing shock that 
is located at the M 87 radio core. In this case, the counter-feature would have to be interpreted 
as the jet before the standing shock. Although the peak of the radio emission is delayed, the 
radio flux started to increase at about the same time as the VHE flare, which indicates that the 
two emission regions are not spatially separated, making this scenario unlikely. B: The VHE 
flare indicates the passage of the disturbance through the standing shock accompanied by a slow 
increase of the radio flux. In this interpretation the first flare related to the injection at the base 
of the jet would have been missed completely in any of the wavelengths. There would still be 
a problem with a non time coincident maximum in the VHE and radio emission. The VHE to 
radio lag could be explained by synchrotron self-absorption (see Sec. 3.), which would seem to 
require a coincidental juxtaposition of the standing shock with just the right radio optical depth 
and subsequent optical depth decline down the expanding jet. An alternative scenario could be 
that the shock-accelerated particles causing the VHE emission cool rapidly until they later emit 
photons dominantly at radio frequencies. 

BLLac is 16 times farther away than M 87, and the jet angle is considerably smaller (9 ~ 
7deg) as compared to M87. The black hole in M87, on the other hand, is ~ 30 times more 
massive than the one in BL Lac. Therefore, our data have a ~ 16 times higher spatial resolutiorj^] 
and provide a more than two orders of magnitude more detailed insight into the jet physics on 
gravitational scales: In case of the M87 observations presented here, lmas corresponds to 
140 i? s , whereas in the case of the BLLac observations, lmas corresponds to ~ 70,000i? s . 
Although the Marscher et al. model makes use of the observed VHE emission in BLLac, 
there is no experimental evidence that would constrain the spatial region of that emission. Our 
observations connect the VHE emission with the radio emission from the nucleus in M 87 and 
therefore contrain the VHE emission region to lie within the collimation region of the jet, at 
maximum a few hundred R s away from the black hole. 
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